D mesons and charmonium states in hot pion matter 
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We calculate the in-medium D meson self-energies in a hot pion gas induced by resonance in- 
teractions with pions. The appropriate resonances in the s, p and d waves of the D meson-pion 
pair are represented by low lying scalar, vector and tensor D* mesons. At temperatures around 200 
MeV the D-meson mass drops by 30 MeV and the scattering width grows up to 60 MeV. Similar 
medium effects are found for the D* vector mesons. This opens and/or enhances the decay and/or 
dissociation channels of the charmonium states \ c and J/9 to DD, D*D, DD* , D* D* pairs 
in pion matter. 
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Since the conjecture of the dissociation of charmonium 
states J/^, . . . in a quark-gluon plasma (QGP) due 
to color screening Q J/$ suppression was considered as 
one of the key signals for the creation of a QGP phase 
in heavy ion reactions. An anomalous J/\I / suppression 
has possibly been seen in central Pb + Pb reactions at 
CERN Super Proton Synchrotron (SPS) 2]. However, 
the interpretation of data is strongly influenced by the 
final state interactions of charmonium states in a hadron 
gas. The dissociation of J/^'s due to pion scattering 
J/* + 7T -> DD, D*D, DD*,D*D* has extensively been 
discussed within the comover scenario 0, for a review 
see 0. The J/^f yield at Relativistic Heavy Ion Col- 
lider (RHIC) will be measured by the PHENIX Collab- 
oration 5]. On the other hand, recent lattice QCD cal- 
culations indicate that the J/9 might exists as a bound 
state even at temperatures above the critical temperature 
for a QGP phase transition 

Direct in-medium modifications of charmonium states, 
i.e. mass shifts and elastic collisional widths are expected 
to be small compared to those of the D mesons Q • Indi- 
rect modifications of charmonium state can, however, ap- 
pear due to medium modifications of the D-meson states. 
With the lowering of D masses the inelastic collisional 
width of charmonium states (due to dissociation to open 
charm) or, in some cases, even the ordinary width (due to 
on-shell decays to open charm) is increased or is appeared 
as a result of the lowering of corresponding thresholds 
These facts are important for both, the production 
and the survival probability of J/9 mesons. About 20%, 
14% and 8% of J/9's are produced in pN interaction 
at 300 GeV/c via the decays of the excited charmonium 
states Xi) X2 and ^'j respectively ||. A dissociation of 
these states in the medium can approximately reduce the 
J/^-production rate by a factor of two. Moreover, the 
shift of the thresholds for the D- meson channels in J/ 9- 
7T collisions increases the J/^-dissociation cross section 
and reduces the J/vP-survival probability. 

D mesons have an obvious analogy to K mesons re- 
placing strange by charm quarks. Medium modifications 



of D mesons and charmonium states have been consid- 
ered by several authors. However, similar as for the kaons 
the investigation of D-meson mass shifts was mainly re- 
stricted to the case of nuclear matter 0, [n], 0, Il3| . 
In ultra-relativistic heavy-ion collisions at CERN SPS 
or RHIC energies the hadronic environment is, how- 
ever, baryon dilute and meson rich EH- Thus, in the 
present work we investigate in-medium properties of D 
and D* mesons and the resulting consequences for char- 
monia in hot pion matter. Previous investigations at 
finite temperature concentrated to large extent on the 
determination of charmonium dissociation cross sections 
0, II El El El El El which are, unfortunately, bur- 
dened by large model uncertainties. 

The strangeness sector has been investigated theoreti- 
cally in [2(j (see also references therein) and it has been 
found that kaons acquire a negative mass shift and a sub- 
stantial broadenin g w hich could explain the enhanced 
d> — > fj, + (i~ decay |21| relative to the 4> — > KK decay 
[22| possibly seen in heavy ion collisions at CERN SPS 
(see although Refs.[23ll24| on possible explanation of the 
effect) and the <f> — > e + e~ to cj> — > KK enhancement 
possibly seen in heavy ion collisions at RHIC 25] . 

The self-energy £ of a D meson in a pion gas is to 
leading order in pion density determined by the invariant 
£>7r-forward scattering amplitude A 
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Scalar n S7r and vector pion densities 
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Here = —\x v - is the n + chemical potential, n^o = 0. 
In isotopically symmetric pion matter which is expected 
to be formed in ultra-relativistic heavy-ion collisions at 
RHIC energies the chemical potentials vanish (n 7T ± = 



0). Hence all pions (negative, neutral and positive) are 
equally distributed. 

In an isotopically symmetric pion gas D + and D° ob- 
tain identical self-energies due to isospin symmetry and 
these are further equal to the D~ and D° self-energies 
due to charge conjugation symmetry. Thus we can re- 
strict ourselves in the following to the discussion of the 
D + self-energy. In this context it should be noticed that 
D mesons show a close analogy to kaons. In both cases 
a mass splitting between K and K, respectively D and 
D occurs at finite nuclear density |12| whereas in isotopi- 
cally symmetric pion matter kaons and antikaons obtain 
equal in-medium modifications |'2C|. 

The interaction of D mesons with pions or more gener- 
ally, the interaction of charmed heavy-light pseudoscalar 
and vector mesons with light pseudoscalar mesons has 
been already discussed in the literature. A local 4- 
particle interaction motivated by chiral symmetry and 
by the similarity of kaons and D mesons was introduced 
in Ref. j2||. This tree level interaction has further been 
iterated for s-wave scattering [2(| and the poles of the 
amplitudes in specific channels were identified with re- 
cently observed scalar and axial D* mesons p?! l28j . The 
tree level £) + 7r-scattering amplitude reads [2(j 
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where F « 93 MeV is the pion decay constant. This 
result can be presented by isospin 1/2 and isospin 3/2 
amplitudes 
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The first amplitude Ai/ 2 is positive and attractive, the 
second amplitude ^3/2 is negative and repulsive. Iterat- 
ing the A1/2 amplitude to all orders in ladder approxi- 
mation leads to a large value for the T matrix (at least 
in s waves) and a pole corresponding to a Dir resonance. 
The A 3 / 2 amplitude after iteration decreases |2(j. The 
tree level amplitudes contain only s and p waves. Not 
discussing the reliability of such an approach in what fol- 
lows we will not rely on these tree level amplitudes but 
instead will take only the observed s- , p- and d- wave 
resonances in the Dtt system into account. For orbitally 
excited mesons we take the experimental data from the 
BELLE Collaboration [2^ (see Fig.l) where for the first 
time all four orbitally excited D-meson states have been 
observed simultaneously. The information on the very 
narrow vector mesons D* (p-wave resonances close to 
threshold) is taken from Particle Data Group[30j. 

For resonance scattering we include only the isospin 
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FIG. 1: (Color online) Levels of D-meson excitations from 
Ref. El. 



TABLE I: Excited D-meson states which are taken into ac- 
count as resonances in Dn system. 



Resonance mass (MeV) width (MeV) 
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and the sum of these amplitudes that enters the ther- 
mal average is equal to Ai/ 2 - For the forward resonance 
amplitude Ai/ 2 we use the relativistic Breit-Wigner form 
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where j = 0, 1, 2 corresponds to the s—,p— and d— wave 
resonances Dq, D* and D 2 in the Dtt system with masses 
Mj , partial and total widths and L' ot , respectively; 
ji = J2 = are the spins of the D and the pion, re- 
spectively, and k is the cm. momentum. The energy 
dependence of the widths is regulated by the specified 
partial wave 
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where subscript '0' refers to on mass shell decay widths. 
When the same amplitude is evaluated for the j% = 1 D* 
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meson, the sum runs over the resonances D[, D±, D\ 
and the ground state D meson. For the branching of 
the tensor resonance we take the world average r(£>2 ~~ * 
Dn) : r(£)| -> D*tt) w 2 : 1. Now the D, respectively 
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FIG. 2: (Color online) The scattering width F = -ImE/M 
and real part of self-energy ReY,/2M as a function of tem- 
perature T are shown for for D mesons (top) and D* mesons 
(bottom) at rest in pion matter . 

the D* self-energy in a pion gas at rest is obtained 
by integration over the pion distribution which yields the 
corresponding modifications of the meson properties: the 
scattering width T and mass shift i£e£/2M. The contri- 
bution of the narrow D* resonance with a width of about 
100 KeV to the D meson medium modifications are of the 
order of the D* width and can be neglected (if one takes 
the in-medium modification of the D* into account, this 
is no more the case as will be seen later on). The contri- 
butions of the scalar and tensor D* mesons are, however, 
sizable and shown in Fig. 2 as a function of the temper- 
ature. Medium modifications of the vector D* mesons 
which are of similar magnitude arise from the coupling 
to axial and tensor mesons, see also Fig. 2. 

For further applications to charmonium states the off- 
shell properties of D and D* mesons, namely the spectral 
functions will be important. Spectral function is defined 
by the off-shell self-energy E(m 2 , |p|) as follows 

i i i i\ 1 — ImS 
P [m , p = -j- — — 2 . 2 j- . 2 ■ 8 



medium with temperature T = 200 MeV the spectral 
functions are shown on Fig. 3 as the functions of the in- 
variant mass. At other momenta typical for the decays 
of the charmonium states (see below) they don't change 
significantly because the velocities of D and D* mesons 
are small compared to the thermal velocities of pions. 
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FIG. 3: (Color online) The in-medium spectral functions of D 
and D* mesons at rest in the hot pion gas with temperature 
T = 200 MeV. The dashed and solid lines for the in-medium 
D,D*-meson spectral functions correspond to the lowest or- 
der calculation and to the first order self-consistent iteration 
of the in-medium spectral functions. 

For D, D* mesons we calculated the spectral functions 
both, with vacuum (dashed lines) and in-medium (solid 
lines) spectral functions of D*,D mesons. In the latter 
case which is the first order self-consistent iteration, the 
vacuum propagators of D and D* in (6) are replaced 
by the corresponding in-medium propagators in order to 
obtain the iterated in-medium results 



S - Ml {D ^ .3 ~ Ml {Dt) - S D(Z 3-) ' 



In medium spectral function depend not only on the in- 
variant mass squared of the particle but also on its mo- 
mentum. For D and D* mesons at rest in the hot pion 



The medium modification of the D meson has thereby 
a feedback on the D* since it modifies the correspond- 
ing decay channel. The difference between the lowest 
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order and the first iteration is significant at small invari- 
ant masses of the D, D* mesons which is important for 
the charmonium decays. 
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FIG. 4: (Color online) In-medium width of J/\l/, Xi> X2 and 
mesons. 



With this information, one is able to evaluate the disso- 
ciation widths of the charmonium states J/*, Xi, X2, 
in the pion medium. Let us take into account the fact 
that in the medium D, D* mesons are not particles 
but resonances with finite collisional widths. Then the 
J/<1> -> DD, DD*, D*D, D*D* decay can occur sub- 
threshold and thus account for the reactions J/^tt — > 
DD, DD*, D*D, D*D* which are open in pion matter. 
The width T( J/* — > DD) can be easily found 



r(j/* -> dd) 
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xp D (rnl,p)p B (rnj,p)dmldmj , 

where p is cm. momentum for the decay of J/ty meson to 
D and D* mesons with masses mi and ma, respectively 
and gj/^DD l& the J/^DD coupling constant. Taking 
the value of gj/^oD = ^.8 from ref. Jg] one obtains a 
dissociation width of r(J/# -> DD) = 0.54 MeV at a 
temperature of T = 200 MeV. The other decay channels 
(DD* D*D, D*D*) can be treated analogously. To do 
so, we use the coupling constants 



gj/tD'D* — 9j/*hdd> 



gj/i>DD* 



9j/9DD 
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in combination with appropriate phenomenological ver- 
tices ^t|. I n total one obtains thus at a temperature of 
T = 200 MeV a J/* collisional width of 1.15 MeV (see 
Fig.4). 

These values can now be compared to the collisional 
width r(J/^7r -► DD, DD*, D*D, D*D*) = 5 
14 MeV. The later one was obtained by multiplying the 
average cross section for the J /iff dissociation through 
pions evaluated at the same temperature and equal to 
< a wJ ^v >« 0.75 -r 2 mb H3 (v being the relative ve- 
locity) with the corresponding pion density. The order of 
magnitude agreement between the J/^f subthreshold de- 
cay width to D mesons and the J/^ir dissociation width 
is not surprising because both effects are of common na- 
ture. In fact, the J/\& dissociation through pion capture 
by one of the D mesons is connected to the D mesons 
collision widths. In the calculation of the J/^> decay rate 
the D mesons collision width appears in the denomina- 
tor of the Breit-Wigner amplitude which, due to the sub- 
threshold character of the J/\& decay, can be expanded 
in powers of the pion density entering the D-meson col- 
lision width. To lowest order the result coincides with 
the J/^f dissociation rate through one-pion capture by 
one of the D mesons. However, the inclusion of realistic 
spectral functions turns out to reduce the J/^f dissoci- 
ation. So, the results of previous calculations [l7j differ 
from the present ones because some contribution (contact 
4-particle interaction) is not considered in our approach 
(see diagram lc in Fig. 1 in .17] ) and because our am- 
plitude is closer to the unitarity constraint. In this con- 
text it should be mentioned that there appear additional 
dissociation mechanisms at the quark level [l8L IT^. l3l| , 
e.g. s-channel box diagrams, which cannot be expressed 
in terms of effective meson degrees of freedom and may 
lead to an additional increase of the width. 

The in-medium width of the excited is estimated 
analogously assuming the same g^ijjf) coupling constant 
as for the J/*. Since the lies only 52 MeV below the 
DD threshold in free space its in-medium width is about 
50 times larger than that of the J /iff. Finally, the in- 
medium widths of the xi, \i states were estimated using 
the phenomenological vertices gMDX^(D fl D + D^D) and 
gMox^D^Dv, g/MoXnud^DduD with g = gj/^nn, 
respectively. As seen from Fig.4 xi, Xi mesons also re- 
ceive a substantial width in pion matter. 

In conclusion, the D meson self-energies in pion mat- 
ter have been determined to leading order in density 
taking thereby resonances in the Dir amplitude into ac- 
count. The resonances in s, p and d waves of D meson- 
pion system were represented by low lying scalar, vec- 
tor and tensor D* mesons which have been observed 
experimentally. This allows to determine the D self- 
energy in a model independent way. At a tempera- 
ture around 200 MeV the D-meson mass is reduced by 
about 30 MeV and the scattering width is about 60 MeV. 
Similar medium modifications were found for D* vec- 
tor mesons. Consequently, the widths of the decay and 
dissociation channels of the charmonium states \ c 
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and J/* to DD, D*D, DD* , D*D* pairs are enhanced 
(from Fj/y ~ 1.15 MeV to IV ~ 62 MeV at T ~ 200 
MeV). As a consequence, feeding of J/ty states from ex- 
cited charmonium states ceases in a hot pion gas which 
characterizes the hadronic final state in high energetic 
heavy ion reactions in good approximation. 

Hence, D mesons modifications in hot pion medium are 
important for the production of J "/W during the fireball 



expansion in heavy ion reactions. The back reactions 
of J/^ formation in D mesons collisions from charmed 
meson rich medium become important for corresponding 
transport simulations. 
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